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Predictions for the charm-production reaction pp — > A^ A+ for energies near the 



threshold are presented. The calculations are performed in a meson-exchange frame- 
work in close analogy to our earlier study on pp — > AA by connecting the two pro- 
cesses via SU(4) symmetry. The obtained A~A+ production cross sections are in 
. the order of 1 to 7 fib, i.e. a factor of around 10 smaller than the corresponding 

cross sections for AA. However, they are 100 to 1000 times larger than predictions 
of other model calculations in the literature. 
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^ ' 1 Introduction 

> 

X 

^ , The study of the production of charmed hadrons in antiproton-proton (pp) 

colhsions is of importance for the understanding of the strong force in the non- 
perturbative regime of QCD. The FAIR project at the GSI laboratory has an 
extensive program aiming at a high-accuracy spectroscopy of charmed hadrons 
and at an investigation of their interactions with ordinary matter [1]. Presently 
httle is known about such interactions and their knowledge is a prerequisite 
for investigating issues like in-medium properties of charmed hadrons, e.g. 
cc-quarkonium dissociation and changes in properties of D mesons due to chi- 
ral symmetry restoration effects on the light quarks composing these mesons. 
Therefore there is an urgent need for theoretical investigations to guide such 
experiments. 

In the present paper we concentrate on the reaction pp A^A+ close to its 
threshold. Providing predictions or simply estimations for this reaction is very 
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challenging. First of all, the lack of any direct empirical information on this 
reaction makes it difficult to constrain model parameters. Second, both long- 
distance and short-distance physics are present in the reaction and this poses 
the question on the appropriate degrees of freedom to be used in order to de- 
scribe it, quarks and gluons or mesons and baryons, or a combination of both. 
To the best of our knowledge, presently there are only four elaborate studies 
that consider this reaction. The most recent publication is by Goritschnig et 
al. [2j, who employ a quark-gluon description based on a factorization hypoth- 
esis of hard and soft processes. This work supersedes an earlier study by that 
group within a quark-diquark picture, where already concrete predictions for 
the A~A+ production cross section were given [3]. In the study by Kaidalov 
and Volkovitsky [4J a non-perturbative quark-gluon string model is used, based 
on secondary Regge pole exchanges including absorptive corrections. On the 
same lines, there is the more recent publication by Titov and Kampfer [S]. 

Our work here builds on the Jiilich meson-baryon model ^6i7j for the reaction 
pp AA. In that model the hyperon-production reaction is considered within 
a coupled-channel approach in which the microscopic strangeness production 
process and the elastic parts of the interactions in the initial (pp) and final 
(AA) states are described by meson exchanges, while annihilation processes are 
accounted for by phenomenological optical potentials. The elastic parts of the 
initial- and final-state interactions (ISI and FSI) are G-parity transforms of an 
one-boson-exchange variant of the Bonn NN potential [8] and of the hyperon- 
nucleon model A of Ref. [S], respectively. The model achieved a reasonably 
good overall description of the wealth of pp AA data collected by the P185 
collaboration at LEAR (CERN) on total and differential cross-sections and 
spin observables [10] - see also the review in Ref. [TT] . 

The extension of the model to the charm sector here follows a strategy sim- 
ilar to our recent work on the DN and DN interactions [T^fTB] , which used 
the original Jiilich meson-exchange model for the KN system [TillTS] and im- 
provements from quark-gluon dynamics at short distances [T6p!7] . Specifically, 
in the present paper we construct an extension of the meson-exchange model 
of Refs. piT] to the reaction pp A^A+ assuming as a working hypothe- 
sis SU(4) symmetry constraints. We examine the sensitivity of the results to 
changes in the elastic and annihilation parts of the initial pp interaction and 
inspect the effects of the final A~A+ interaction and of the form factors that 
enter in the pp A~A+ transition potential described by t-channel D and D* 
meson exchanges. In addition, we also investigate the effect of replacing this 
meson-exchange transition by a charm-production potential derived in a quark 
model. We believe this is important for assessing uncertainties in the model, 
since one could easily raise the question on the validity of a meson-exchange 
description of the transition in view of the large masses of the exchanged 
mesons. 
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In the next Section we discuss the basic ingredients of the model and fix pa- 
rameters by fitting differential and total inclusive pp cross sections. In Section[3] 
we present numerical results for our predictions for differential and total cross 
sections and compare with the results available in the literature. A summary 
of our work is presented in Section HI 



2 The model 



We will start discussing the basic ingredients of the original Jiilich coupled 
channel approach Pf7] that we employ here for the reaction pp —>■ A^rA^. 
The transition amplitude is obtained from the solution of a coupled-channel 
Lippmann-Schwinger equation, 



T(q', q, z) = V{c^, c,, z) + j d'q" V{c{, q", z) Go(q", z) T(q", q, z), (1) 

where z is the initial energy and q' (q) the cm. (center-of-mass) relative 
momentum in the initial (final) state. Here ^ is a 2 x 2 matrix in channel space 
containing the interaction potentials (channel 1 = pp, channel 2 = A~A^) 



y(q',q,^) 



v2Hq',q,^) v^''(q',q,^) 



(2) 



and Go(q, z) is the propagator 
/ 



G'o(q,^) 



l/(z-E«+^e) 

l/{z-E^^^ + ie) 



(3) 



with = E^{p) + E^{p) and E^^'> = E^{Ac) + E^{Ac). The diagonal poten- 
tials V" are given by the sum of an elastic part and an annihilation part. 

Though the Jiilich group has developed rather sophisticated models of the 
NN interaction p^jTO] these potentials cannot be used anymore at such high 
energies. Indeed, already in the study of pp AA the elastic part of the 
NN potential was deduced (via G-parity transform) from a simple, energy- 
independent one-boson-exchange A^A^ potential (OBEPF) and a phenomeno- 
logical spin-, isospin-, and energy-independent optical potential of Gaussian 
form, 

VX''{r) = {Uo + tWo)e-^'/'^i (4) 
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Fig. 1. Differential cross section for elastic pp scattering at piab = 10.1 GeV/c as a 
function of t. The dash-dotted curve corresponds to a calculation where only one 
pion exchange is added to the optical potential (A). The dotted curve results when 
the complete G-parity transformed OBEPF model of Ref. [8j is used for the elastic 
part (D). The dashed and solid curves are obtained by leaving out vector-meson 
exchanges (B) or by reducing the elastic part (except for the pion exchange) to 
10 % (C), respectively. The experimental information is taken from Foley et al. [2l] 
and Berglund et al. [23] . 

was added in order to take into account annihilation. Now, at even much 
higher energies, any A^A^ potential has to be considered as being purely phe- 
nomenological. Still we keep the longest ranged (and model-independent) part 
of the elastic pp interaction, namely one-pion exchange. To it we add again an 
optical potential of the form given in Eq. (jlj) and determine the parameters 
{Uq, Wo, To, cf. Table [I]) by a fit to NN data in the energy range relevant for 
the reaction pp A^A^. Fortunately, there are total cross sections [20|21f22] 
and even differential cross sections [2T|23] around piab = 10 GeV/c, i.e. fairly 
close to the A~A+ threshold which is at 10.162 GeV/c. This pp potential is 
called model A in the following. As can be seen from Table [2] and Fig. [1] the 
integrated cross sections as well as the pp different cross section are fairly well 
reproduced. 

Since one knows from studies on pp —>■ AA that the magnitude of the cross 
sections depends very sensitively on the ISI p[7|24|25] we consider here several 
variants of the NN potential. First, we take the full G-parity transformed 
OBEPF interaction and add again the optical potential (model D). Obviously 
in this case, only the total cross section can be still brought in line with the 
experiment by adjusting the parameters of the optical potential, while all other 
cross sections are strongly overestimated. Still we consider this potential too 
with the intention that it serves as an illustration for the uncertainties in our 
predictions due to the used pp interaction. Realizing that the overestimation 
of the cross sections is primarily caused by the vector-meson exchange (p, u) 
contributions ^26j to OBEPF, we prepare two more potentials where we (i) 
leave out those vector mesons altogether (B) or (ii) reduce the elastic part 
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(except for one-pion exchange) to 10% (C). Within both scenarios a rather 
satisfying description of the NN data around 10 GeV/c can be obtained, cf. 
Tableland Fig.[TJ In particular, not only the slope but even the shoulder in the 
differential cross section is reproduced quantitatively by these two interactions. 

Table 1 

Parameters of the phenomenological optical potential in the pp channel for the four 
different models described in the text. 





A 


B 


C 


D 




-48.0 MeV 


100.5 MeV 


72.9 MeV 


1808 MeV 


Wo 


-531.9 MeV 


-529.5 MeV 


-448.3 MeV 


-1644 MeV 


ro 


0.56 fm 


0.56 fm 


0.58 fm 


0.41 fm 



Table 2 

Integrated pp cross sections at piab = 10.1 GeV/c for the four models considered. 
The experimental value for the charge-exchange channel pp — > nn (cex) is for piab = 
9 GeV/c [22] . 





crtot [nib] 


aei [mb] 


acex [mb] 


experiment 


54.7 ± 0.60 [20J 


14.6 ±3.3 [21J 


0.284 ±0.041 [22] 


A 


54.2 


14.4 


0.20 


B 


54.1 


14.6 


0.23 


C 


54.3 


14.2 


0.42 


D 


54.6 


20.6 


2.82 



The interaction in the final A^A+ system is assumed to be the same as the one 
in A A. Specifically, this means that the elastic part of the interaction is given 
by the isospin-zero a and u exchanges with coupling constants taken from 
the hyperon-nucleon model A of Ref . [2] , while the annihilation part is again 
parameterized by an optical potential which contains, however, spin-orbit and 
tensor components in addition to a central component [6j: 

yKAt^K^t^^^ = [U, + iW, + {Uls + ^WLs) L ■ S 

+ iUt + iWt)aA^-ra-^^-r]e-''"/''o. (5) 

The free parameters in the optical AA potential were determined in Ref. [6] by 
a fit to data on total and differential cross sections, and analyzing power for 
pp — > A A. Clearly, there are no reasons to believe that the A~A+ interaction 
will be the same on a quantitative level. But we expect that at least the bulk 
properties are similar. Specifically, in both cases near threshold the interac- 
tions will be govered by strong annihilation processes. We will investigate the 
role played by the FSI for our results by simply switching it off. We use the 
parameters as given in Table II of Ref. |7j. 

In principle, the parameters of V^t~^^^ can only be determined together with 
those for A~A+ because we consider a coupled-channels problem. However, 
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like in AA the branching ratio pp A^ A^ is so small that the effect of the 
A~A+ channel on the diagonal pp T-matrix can be savely neglected. 

The transition potential from pp to A^A+ is given by t-channel D and D* 
exchanges and their explicit expressions are the same as for K and K* as 
given in Ref. [H]. They are of the generic form 

yPP^KAt^t)^ E 9%A.M^f^, (6) 
M=D,D* ^ "^M 



where QnAcM are coupling constants and -^atAcA/I^) form factors. Under 
the assumption of SU (4) symmetry the coupling constants are the same as in 
the corresponding exchanges in pp AA. However, the cutoff masses in the 
vertex form factors cannot be taken over, because the masses of the exchanged 
particles are now much larger. We use here a monopole form factor with a 
cutoff mass A of 3 GeV, at the NAcD as well as at the NAcD* vertex but we 
will explore the sensitivity of the results to variations of the cutoff mass. 

To examine further the uncertainties, as an alternative we consider here also 
a charm-production potential inspired by quark-gluon dynamics. There is a 
large literature associated with strange-hadron production in pp reactions, the 
best known works are those of Kohno and Weise [23] , Furui and Faessler 
Burkardt and Dillig [2E], Roberts [2S| and Alberg, Henley and Wilets 
a more complete list of references can be found in the review of Ref. [TT] . In 
the present study we adopt the interaction proposed by Kohno and Weise, 
derived in the so-called ^Si mechanism of a constituent quark model. In this 
model the ss pair in the final state is created from an initial uu or dd pair via 
s-channel gluon exchange. After quark degrees-of-freedom are integrated out 
the potential has the form |24j: 

V^-^\r) = frr^Ss,6ro [^^j x exp(-3rV(4(r^))) . (7) 



Here a/rriQ is an effective (quark-gluon) coupling strength, (r^) is the mean 
square radius associated with the quark distribution in the p or A, and S and 
T are the total spin and isospin in the pp system. This simple potential has 
actually a very modern form, namely that of a contact term, though smeared 
out by the quark distribution. The effective coupling strength is practically a 
free parameter that was fixed by a fit to the pp — > AA data [7] . But it depends 
implicitly on the effective gluon propagator, i.e. on the square of the energy 
transfer from initial to final quark pair, cf. Refs. [27ll28f30j . Heuristically this 
energy transfer corresponds roughly to the masses of the produced constituent 
quarks, i.e. mc ~ 2mg, we expect the effective coupling strength a/rriQ for 
charm production to be reduced by the ratio of the constituent quark masses 
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of the strange and the charmed quark squared, {ms/rricY ~ (550 MeV / 1600 
MeV)^ ~ 1/9 as compared to the one for pp — » AA. Note that a different 
suppression factor arises for the ^Pq quark-antiquark annihilation vertex, con- 
sidered by Alberg et al. in their study of pp AA [23] • In this case the 
amphtude scales with l/rriq so that the effective strength of the transition 
potential is reduced by rris/mc ~ 1/3 only when going from strangeness to 
charm production. 



3 Results 

Our predictions for the total reaction cross section for pp ^c^t ^'^^ P^^" 
sented in Fig. [2l The cross section is shown as a function of the excess energy 
e = ^/s — — mx^ so that we can compare it with the one for pp —>■ AA 
at the corresponding e. The curve in Fig. [2] correspond to the pp model C, for 
which we obtained the largest results. 

Obviously, and as expected, the cross section for A^A+ production is smaller 
than the one for AA. But the difference is about one order of magnitude 
only. Indeed, one can get a very rough estimation for the reduction just by 
considering the following: When going from AA to A^A^t within the meson- 
exchange picture the main change occurs in the meson propagators of the 
transition potentials where the masses of the strange mesons {K, K*) are re- 
placed by those of the charmed mesons (Z), D*), niM^ — ^ f^M^ cf. Eq. (|6]). 
All coupling constants remain the same because of the assumed SU(4) sym- 
metry. Thus, the ratio of the transition potentials is then given roughly by 
ypp^kckc lypp^AK ^ rn\jjm\j^ ~ 1/4 so that one expects the cross section 
to be smaller by a factor of around 16. Our explicit calculation is pretty much 
in line with this admittedly rather qualitative estimation. 

Let us now discuss the sensitivity of the results to the various ingredients of our 
model calculation. The most crucial component is certainly the pp interaction 
in the initial state. Without it, and specifically in Born approximation, we 
would get cross sections that are more than a factor 100 larger. This is no 
surprise because the very same gross overestimation in the Born approximation 
occurs for pp —>■ AA [6H71l24f25] but also for other exclusive pp annihilation 
channels like pp — >• tttt, pp — >• KK, etc. [TSfTQ] . if one does not take into 
account the ISI. However, with the ISI included, it is reassuring to see that 
the variations of the predicted pp A~A+ cross section for the different pp 
potentials we have prepared is fairly small. To be concrete, the cross sections 
at e = 100 MeV are 5.8, 6.3, and 7.3 fib for the interactions A, B, and C, 
respectively. Only the result based on the full G-parity transformed OBEPF 
potential differs more significantly. Here we get a cross section that is with 0.8 
nb a factor 10 smaller than the other values. The suppression for this pp model 
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Fig. 2. Total reaction cross sections for pp A A and pp — > A^A+ as a function of 
the excess energy e. The results for pp — > AA (upper curves) are taken from our work 
[7]. The solid curves are results for the meson-exchange transition potential while 
the dashed curves correspond to quark- gluon dynamics. The pp A~A^ results 
are obtained with the pp interaction C. 
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Fig. 3. Differential cross section, singlet fraction, depolarization Dnn and spin-cor- 
relation parameter Cyy for pp — > A~A+ at piah = 10.343 GeV/c (e = 40 MeV) as 
a function of the the cm. scattering angle. The curves correspond to different ISI 
as indicated in the figure. The meson-exchange potential is used for the transition 
interaction. 

is presumably due to the fact that it yields a elastic cross section that is 
much larger than the ones of the other potentials considered, and actually in 
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disagreement with experimental information. Still, the variation of one order 
of magnitude may be considered as a realistic estimation of the uncertainty 
due to the ISI. 

Once the ISI is included our results turned out to be rather insensitive to the 
final A~A^ interaction. Even when we leave it out altogether the cross sections 
do not change (decrease) by more than 10-15 %. We also considered variations 
of the cutoff mass at the NA^D and NAcD* vertices. For the results discussed 
above a value of A = 3 GeV has been used. When reducing this value to 
2.5 GeV the cross sections at e = 100 MeV drop by roughly a factor 3. Of 
course, employing even smaller cutoff masses would further decrease the cross 
sections. However, since the exchanged mesons have a mass of around 1.9 to 
2 GeV we consider values below 2.5 GeV as being not really realistic. 

We display here also the results based on an adaption of the ^Si quark-gluon 
transition mechanism of Ref. [23]. We scale the effective coupling strength 
a/rriQ = 0.25 fm^, fixed in our study of pp — > AA [7], with [rric/msY using the 
constituent quark masses = 550 MeV and rric = 1600 MeV, i.e. the same 
values as employed in our previous works in Refs. [T2II3]. As expected, we 
obtain cross sections that are of the same magnitude as those predicted in the 
meson-exchange picture though roughly a factor three smaller, cf. the dashed 
line in Fig.[2l In principle the factor (r^) in Eq. ([7]) also changes when replacing 
nis by rric, since it gives the size of the quark distribution of the overlapping 
hadrons - (r^) decreases as is replaced by rric. We have not changed this 
factor in the present work, since we expect this change to be much smaller 
than the change in the square of the energy transfer m^, since hadron sizes do 
not scale as simple powers of quark masses and so the overlap between hadron 
wave functions is not expected to change drastically from A to Ac - for explicit 
numbers for the size parameters of harmonic oscillator wave functions for N 
and Ac see Ref. [31] • The size parameter we use is (r^)^/^ = 0.55 fm. 

Predictions for differential cross sections, the singlet fraction Fg, the depolar- 
ization Dnn and the spin-correlation parameter Cyy (cf. Ref. [7] for definitions) 
at e = 40 MeV are presented in Fig. [3] for all four ISI considered. Note that the 
singlet fraction is defined by = — ■ (^2)) = ^(1 + Cxx — Cyy + C^z) [6J. 
It is zero if the A~A+ pair is produced purely in a triplet state. Indeed in case 
of pp — i> AA it was found experimentally that the singlet fraction is close to 
zero [lOj and this feature was reproduced by meson-exchange potentials [6] as 
well as by transition interactions based on quark-gluon dynamics [7ll24f25j . In 
the former it is due to a strong tensor force generated by the combined K + K* 
exchange which leads to a dominance of transitions of the form L^^ = Lpp — 2 
where L refers to the orbital angular momentum. As can be seen from Fig. [3l 
the very same feature is predicted also for the reaction pp A~A^. The dif- 
ferential cross section exhibits a peak in forward direction, similar to what has 
been also observed for pp — > AA piTpiI] . 
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p„,(GeV/c) 

Fig. 4. Total reaction cross sections for pp —>■ A~ A+ as a function of the excess energy 
piab- The dark (red) shaded band (blue grid) is the prediction of our meson-exchange 
(quark-gluon) transition potential. The dotted curve is the result from Ref. [1] while 
the dash-dotted curve and the corresponding (green) band is from Ref. [2]. 

Let us now compare our predictions with those by other groups. This is done 
in Fig. m Goritschnig et ah [2j as well as Kaidalov and Volkovitsky [1] have 
presented explicit results in their publications and we reproduce them in Fig. H] 
to facilitate a comparison. Our results are shown as dark (red) shaded band 
(and grid) in order to reflect the variation of the predictions when the four 
different ISI are used. It is remarkable that our results differ drastically from 
those of the preceeding works. Specifically, our cross sections are a factor 1000 
larger than those given by Goritschnig et al. and they are still about 100 times 
larger than the ones by Kaidalov and Volkovitsky. Thus, even when considering 
the variation of about a factor ten due to the ISI that we see in our results 
and the uncertainties due to the unconstrained FSI and form factors in the 
transition potential that amount to roughly a factor three, we are faced with 
an impressive qualitative difference. 



4 Summary 

In this paper we presented predictions for the charm-production reaction pp —>■ 
A~A^. The calculations were performed in the meson-exchange framework in 
close analogy to our earlier study on pp —>■ AA by connecting the two processes 
via SU(4) symmetry. The interaction in the inital pp interaction, which plays 
a crucial role for the quantiative predictions and which is now needed at a 
much higher energy, is re-adjusted so that available pp scattering data in the 
relevant energy range are reproduced. 

The obtained A^A+ production cross sections are in the order of 1 to 7 fj,b. 
Thus, they are just about a factor 10 smaller than the corresponding cross 
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sections for AA. The reduction is in line with a naive estimation based on 
the mass ratio of strange {K, K*) versus charmed {D, D*) mesons. The 
exchange of those mesons governs the range of the forces that are responsible 
for producing a A or a A^,, respectively, in the meson-exchange picture. 

Interestingly, other model calculations in the literature predict A~A+ produc- 
tion cross sections that are 100 to 1000 times smaller than what we found. 
Since the approaches used in those calculations are rather different from ours 
it is practically impossible to say where this drastic difference in the predicted 
cross sections comes from. Hence, a discrimination between those scenarios 
appears to be possible only when the reaction pp — > A~A+ will be examined 
in a future experiment at FAIR. Of course, if our predictions turn out to be 
closer to reality than the others, it will be much easier to perform a pertinent 
experiment at FAIR. Moreover, Ac's and Ac's can then be produced more co- 
piously so that also any other experiment that requires charmed baryons as a 
probe will be much more promising. 

The model presented here leaves room for improvements in several directions. 
In particular, it would be important to extend the simple picture of s-channel 
qq annihilation followed by cc creation for the transition amplitude in the 
quark model as used here. An interesting possibility is the one pursued in 
Ref. [29] for pp — >■ AA where interacting many-quark intermediate states are 
considered in the process. 
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